Thermal-radiation-induced nonequilibrium carriers in an intrinsic graphene 
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We examine an intrinsic graphene connected to the phonon thermostat at temperature T under 
irradiation of thermal photons with temperature T r , other than T. The distribution of nonequi- 
librium electron-hole pairs was obtained for the cases of low and high concentration of carriers. 
For the case when the interparticle scattering is unessential, the distribution function is determined 
by the interplay of intraband relaxation of energy due to acoustic phonons and interband radia- 
tive transitions caused by the thermal radiation. When the Coulomb scattering dominates, then 
the quasi-equilibrium distribution with effective temperature and non-equilibrium concentration, 
determined through balance equations, is realized. Due to the effect of thermal radiation with tem- 
perature T r 7^ T concentration and conductivity of carriers in graphene modify essentially. It is 
demonstrated, that at T r > T the negative interband absorption, caused by the inversion of carriers 
distribution, can occur, i.e. graphene can be unstable under thermal irradiation. 

PACS numbers: 73.50.Fq, 73.63.-b, 81.05.Uv 



Different kinetic phenomena caused by carriers local- 
ized near the band cross-point of graphene, including 
dc (magneto)transport and optical properties, have been 
studied intensively within recent years, see reviews [l| 
and last references in 0, H|. The main attention was 
paid to examination of linear response of the carriers in 
the phonon thermostat at the temperature T. Because 
of the weak carrier interaction with acoustic phonons [4] 
different external factors can easily disturb the equilib- 
rium of electron-hole system, and the linear responce be- 
haviour realization needs accurate control. In particular, 
when the sample is not isolated from external thermal 
radiation with temperature T r ^ T, the carriers inter- 
action with additional thermostat of thermal photons is 
essential. This interaction can be effective enough, be- 
cause the interband transitions are determined by the 
velocity vw — 10 8 cm/s characterizing the linear spec- 
trum of carriers (the neutrino-like states near the band- 
crossing point are described by the Weyl- Wallace model 
PJ). Therefore graphene is very sensitive for thermal 
irradiation: the concentration and conductivity of carri- 
ers modify essentially (particularly, the photoconductiv- 
ity induced by thermal irradiation occurs, compare with 
|6|, where the case of the interband pumping was dis- 
cussed). 

In this paper, the results for the distribution of non- 
equilibrium carriers in the intrinsic graphene, interact- 
ing with phonon and proton thermostats with different 
temperatures, are presented. This distribution is ob- 
tained from the kinetic equation, taking into consider- 
ation the quasi-elastic energy relaxation due to acous- 
tic phonons, and generation-recombination processes due 
to interband transitions, caused by thermal irradiation 
(the corresponding collision integrals were obtained in 
Q). Under the high concentrations it is also necessary 
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to take into consideration the Coulomb scattering, which 
does not cause the interband transitions, see p}. More- 
over, the scattering by static disorder should be taken 
into account as the main mechanism of momentum re- 
laxation Q. In the low temperatures range, where the 
carrier concentration is not high, the Coulomb scatter- 
ing is unessential and the distribution function differs es- 
sentially from equilibrium one due to interplay between 
acoustic scattering and radiative transitions. At high 
temperatures, when the Coulomb scattering dominates, 
the quasi-equilibrium distribution of carriers, with effec- 
tive temperature and non-equilibrium concentration, is 
imposed. We also calculate the dependences of concen- 
tration and conductivity on T and T r under the scatter- 
ing by the short-range static disorder. 

In the intrinsic graphene with the symmetrical c- and 
w-bands, and with similar scattering in these bands, the 
distributions of electrons and holes are identical; they 
are described by distribution function f p . This function 
is governed by the quasi-classic kinetic equation [6j]: 

J&LA(f\p) + Jw(f\p) + hc{f\p) = 0. (1) 

Here the collision integrals Jqla, Jsr and Jqc describe 
the relaxation of carriers caused by the phonon (LA) 
and photon (i?) thermostats and the carrier-carrier scat- 
tering (C), respectively. The solutions of Eq (1) have 
been obtained below for the two cases: (a) low concen- 
trations, when Jbc can be neglected, and (b) high con- 
centrations, when Jfic imposes the quasi-equilibrium dis- 
tribution with parameters determined from the equations 
of the balance of concentration and energy. After sum- 
mation of Eq. (1) over p-plane with the weights 1 and p 
we get these balance equations in the form [9(: 

-7^E J ^IP)= ' ( 2 ) 
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where the contribution of acoustic scattering is omitted 
from the equation of the concentration balance, because 
the interband transitions are forbidden due to the in- 
equality s <C vw ■ The inter-carrier scattering does not 
change the concentration and energy. [7[ 

We start with the examination of low temperature 
case, when the carriers concentration is small and Jec 
in Eq.(l) can be neglected. With the use of the collision 
integrals Jrla and Jrr, presented in with temper- 
atures T and T r correspondingly, we get the non-linear 
equation of the second order for the distribution function 
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Here we introduce the characteristic momenta px — 

T/v w and p r = T r /v w ; N 2p / Pr = [exp(2p/p r ) - 

is the Plank function. The rates of quasi-elastic relax- 



ation at acoustic phonons 

B(r) 
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p — (s/vw) v ac p/h, and 
radiative transitions, v p y ' ' — v r p/h, have been expressed 
through the sound velocity, s, and through the character- 
istic velocities, v ac , and v r . They separate the momen- 
tum dependence of the relaxation rates, which is propor- 
tional to the density of states. According to Q, where 
the temperature dependence of mobility have been ex- 
amined, v ac — 1-35 • 10 cm/s for the helium temperature 
(moreover v ac oc T), and the quasi-elastic character of 
scattering is determined by the small parameter s/vw- 
For the case of graphene, placed between the SiC>2 sub- 
strate and cover layer, we get v r ~41.6 cm/s, see @. 

The boundary conditions in Eq.(4) are imposed both 
by the demand of the finite stream along the energy axis, 
when p 4 (dfp /dp + f p ) p ^ 00 < const, and by the concen- 
tration balance equation (2). Because the acoustic and 
radiative contributions become zero under the equilib- 
rium contributions with temperatures T and T r corre- 
spondingly, and these two contributions trend to 1/2 at 
small p, we get f p ^o — 1/2- 03 This demand can be 
used as the boundary condition for the numerical solu- 
tion of Eq.(4). Note, that the validity of Eq. (2) for the 
obtained distribution should be examined. This solution 
was carried out below with the use of the finite difference 
method and the iterations over non-linear contributions 
in Eq.(4), see [ll|. 

In Fig. la-c the obtained distribution functions versus 
energy vyyp for different T and T r are presented. Due 
to the smallness of acoustic contribution at p — ► (see 
above), f p is close to equilibrium distribution with the 
temperature T r . If T r > T the distribution increases at 
high energies up to the range, where Jla is dominant. 
With the further increase of vwP the distribution de- 
creases rapidly on the scale of energies T, so that the 
peak of distribution is formed at the energies ~ T r ; this 
peak causes the carriers concentration of the order of 
equilibrium value at the temperature T r . In this case the 
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FIG. 1: (Color online) Distribution functions governed by Eq. 

(4) for the cases: T=4.2 K (a), T=77 K (b), and T=150 K 
(c) at T r =4.2 K (1), 20 K (2), 77 K (3), 150 K (4), and 250 K 

(5) . Dotted curves are the equilibrium thermal distributions 
at T = T r . 



condition f max > 1/2 is realised, i.e. the inverse distribu- 
tion of the carriers occupation takes place for the energies 
close to the maximum of distribution (see below) . On the 
contrary, if T r < T, the equilibrium contribution of the 
slow carriers is replaced by the rapidly decreasing part in 
the range of vwP < T. This distribution determines the 
small concentration of carriers in comparison with the 
equilibrium value at the temperature T. 

Later we shall examine the case of high temperatures 
(and concentrations), when the Coulomb scattering dom- 
inates, imposing the quasi-equilibrium distribution 



f p = {evp[(v w p - fi)/T c ] + iy 



(5) 



The effective temperature of carriers T c and the chemi- 
cal potential n in this distribution are obtained from the 
balance equations (2) and (3). After introducing the di- 
mensionless momentum x = vwp/T Cl we get the equation 
of the concentration balance 
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which imposes the relation between a, and T c /T r , while 
the phonon temperature is omitted out of this equation. 
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FIG. 2: (Color online) Effective temperature T c (a) and maxi- 
mal distribution / p =o (b) versus T r for different temperatures 
T =77 K (I), f50 K (2), 250 K (3), 350 K (4) and 450 K (5). 

In these variables the equation of the energy balance can 
be presented as: 



1 - 2/. 



dxX I e 2xTJT r 



T r — T 



-7" 



T 



dxx A e x -^~f x = 0, 



(7) 



where 7 = (s/vw) 2 v ac /v r oc T determines the relative 
contribution of the phonon and photon thermostats. 

The solution of the transcendental equations (6) and 
(7) gives the distribution (5), dependent on T and T r , 
which can be characterised by the effective tempera- 
ture T c , and the maximum value of the function f p= o — 
[exp(— fJ./T c ) + These values are presented in Fig. 2 

for the concentrations, greater than 3.5 TO 10 cm -2 . Note, 
that the temperature T c differs from T unessentially un- 



der the great change of T r , meanwhile /i and f p modify 
essentially for the temperature range up to 450 K (the 
high-temperature measurements of graphene were car- 
ried out recently [Hj]). The essential difference of the 
case under examination from the low temperature solu- 
tion, presented in Fig.l, is that the value f p —o is being 
fixed (and equal to 1/2), while the value f p= o decreases 
with the increase of T, and increases with the increase of 
T r : note, that f p = =l/2 when T = T r . 
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FIG. 3: (Color online) Carrier concentrations given by Eq. (8) 
versus temperature T r for different T. Solutions of Eq. (4) 
are plotted in the range n < 5 ■ 10 10 cm~ 2 (dashed curves). 
Solutions of the balance equations (6) and (7) are plotted in 
the range n > 3.5 • 10 10 cm" 2 (solid curves). The equilibrium 
concentration (if T — T r ) is shown as dotted curve. 

As one can see from Figs. la,b or 2b, the distribu- 
tion f p or f p can be greater than 1/2 in a certain en- 
ergy range or at low energies. Since the inversion of 
electron-hole pairs occupation, the regime of the neg- 
ative interband absorption can occur because the real 
part of dynamic conductivity is given by expression Q 
Rea^j = (e 2 /4fi,)(l — 2/ Plu ). In the range of parameters 
under examination the intrinsic graphene tends to be un- 
stable, if an additional adsorption is weak enough. 

The distributions, presented in Figs. 1 and 2, are valid 
for the cases of low and high concentrations, correspond- 
ingly. For the calculation of non-equilibrium concentra- 
tion n, dependent on T, and T r , we use the standard 
expression 



00 



dxxf x , (8) 



where the right equality was written under the substitu- 
tion of Eq.(5) into the standard formula. In Fig. 3 we 
plot n versus T and T r for the low- and high temper- 
ature regions, when the Coulomb scattering can cither 
be neglected, or it dominates. At T — T r these curves 
intersect with equilibrium concentration oc T 2 : because 
the concentration is controlled by thermal irradiation, n 
is smaller (or greater) then equilibrium concentration at 
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T < T r (or T > T r ). Note, that the dependences ob- 
tained from equation (4) , and from the equations of bal- 
ance (6), (7) correspond well in the range of intermediate 
concentrations, 3.5 -j- 5 • 10 10 cm~ 2 . 
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FIG. 4: (Color online) (a) Normalized conductivity, a/a , 
versus T r for different T at l c =20 nm. (b) The same for 
T=150 K at different correlation lengths, l c . 

The modifications of the distribution of non- 
equilibrium carriers in the intrinsic graphene, described 
above, lead to modification of conductivity, a. For the 



scattering of momenta due to static disorder with corre- 
lation length l c we calculate the conductivity, using the 
formula [9j 
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where ^(z) — exp(— z 2 )I\{z 2 ) j z 1 is written through the 
first order Bessel function of the imagined argument, 
Ii(z), and a Q is the conductivity in the case of short- 
range disorder scattering, when l c = 0. For the case of 
short-range scattering, pl c /h «C 1 (p is the characteris- 
tic momentum of non-equilibrium carriers) the conduc- 
tivity is written through the low-temperature distribu- 
tion: cr ~ 2<7 / p= o. For the low temperature range, when 
f P =o = 1/2, the conductivity depends weakly on T, and 
T r . On the contrary, in the high temperature range the 
temperature dependences of a on T and T r are essential, 
see Fig. 4a, and Fig. 2b, where / P =o is plotted. With 
the increase of l c the second term in Eq. (9) becomes es- 
sential and the thermal dependences o j o Q become much 
stronger, see Fig. 4b. 

Next, we list the assumptions used. The main restric- 
tion is the examination of the limit cases either of no 
inter-carrier collisions, or of their domination only. De- 
spite some results correspond well in the intermediate 
range of concentrations and temperatures (a disagree- 
ment of ct/ct , versus T r takes place at l c <20 nm only, 
see Fig. 4b), the accurate analysis of the intermediate 
range is beyond of the scope of this paper. Other as- 
sumptions, such as the models of energy spectrum, or 
the scattering mechanisms, as well as the simplifications 
of distribution function used, are rather standard for the 
calculations of the transport phenomena. 

In closing, the consideration performed demonstrates 
the essential effect of thermal irradiation on the intrin- 
sic graphene properties, therefore the transport measure- 
ments should be carried under the control of the condi- 
tion T = T r . Within the study of the device applications 
we should also consider the possible difference between 
the temperature of phonons and the temperature of the 
external thermal radiation. 
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